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Abstract 


We construct a class of representations of the Heisenberg algebra in terms of 
the complex shift operators subject to the proper continuous limit imposed by 
the correspondence principle. We find a suitable Hilbert space formulation of 
our construction for two types of shifts: (1) real shifts, (2) purely imaginary 
shifts. The representations involving imaginary shifts are free of spectrum 
doubling. We determine the corresponding coordinate and momentum oper¬ 
ators satisfying the canonical commutation relations. The eigenvalues of the 
coordinate operator are in both cases discrete. 
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1. Introduction. The aim of this paper is to find a possibly wide class of difference 
operator representations for the coordinate and momentum operators leading to a general¬ 
ization of the standard quantum mechanics (QM) in a way compatible with the discreteness 
of space-time. It is important to stress that at this stage we are at the level of kinematics. 
Thus, we are looking for proper dehnitions of the coordinate and momentum operators, 
their domains, the ’’classical” limit” (to recover the standard QM in the continuous limit), 
the Hilbert space etc. We assume that the Heisenberg commutation relations (CR) hold 
on some dense set of states. A deeper understanding of, seemingly, exotic realizations of 
CR for which one of the operators is bounded might be needed to advance this program. 
In other words, the lack of such understanding was perhaps responsible for the relatively 
unsuccessful, albeit very interesting, results of the early attempts in this direction 0-0. 

Our research is motivated by an old problem of the nature of the space-time geometry 
below the Planck scale. There have been many speculations and indications that at this level 
new discrete structures are likely to emerge (see the arguments put forward in string theory 
0], in various approaches to quantum gravity |]6|-0pT|, in’t Hooft’s approach [|T0| etc.) and 
that, as a consequence, some ’’unusual” representations of the Heisenberg algebra maybe 
appropriate |TT[. Finally, discrete models are widely used because of technical reasons. 


perhaps the most obvious examples beeing lattice gauge theories 
approach to general relativity [O] to name only a few. 


and Regge calculus 


The hrst step to realize this program has been taken in |]T0| , where representations of the 
Heisenberg CR in terms of real shift operators have been analyzed. We have discussed the 
Bargmann-Fock (annihilation-creation operators) as well as the Schrodinger (position and 
momentum operators) representations. In particular, we have found that our coordinate op¬ 
erator has a discrete spectrum. In this paper we generalize our previous results by including 
the complex shift operators as well. In particular, we find that the operators discussed in 
ITsHFTI are included in our scheme. Those operators correspond to purely imaginary shifts 
and in the last two references they are used in a physically interesting context of the black 
hole’s energy quantization. 

We work directly with the coordinate (Schrodinger) representation, as opposed to the 
annihilation-creation type, for the former seems to be physically more fundamental; the 
sentiment captured so nicely in the following quote “In physics the only observations we 
must consider are position observations...' 


2. Momentum and coordinate operators. We start with the following Ansatz for the 
generalized momentum operator P and the corresponding “discrete” derivative D 


N 


P = —ihD = —ih -— ttfc E 

where the complex shift operator is dehned as 

fi^) = f{x + kAx) , 


k 


( 1 ) 


( 2 ) 


and the shift Ax can be complex. We keep the notation the same as in [14| with the 


exception of the overall 1/Ax factor in the formula (0) that dehnes the constants a^’s. With 
this dehnition our formal equations for become identical to those of [Q. The coordinate 
operator X is assumed to have the form 
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( 3 ) 


X = 




xEL + E. 


Ax 


X 


with unspecified constants f3k- 

We would like to have the standard coordinate {x) and momentum {p = —ihd/dx) in 
the continuum limit. Hence, we impose the classical limit conditions 


lim X = X , 

Aa:—>0 

In addition, we demand the symmetry: 


lim P = p 

Ax^O 


( 4 ) 


a_k = oik (5) 

and that D be the best fit to d/dx, i.e. we assume the operator D to be optimal in the 


sense of [Q. This determines the coefficients ak in a unique way, independent of Ax 

{N\f 


^ + k)\{N - k)\ ' 

Finally, the operator X will be determined from the Heisenberg CR 

[X, P] = ihl . 


( 6 ) 


( 7 ) 


At the formal level imply the same solutions for (5^8 as in the case of the real 
shifts |n]. However, to determine Hermicity we have to specify the scalar product and the 
corresponding Hilbert space. This is an especially pressing issue for the imaginary shifts 
where the standard setup involving functions in L^(IR) is not suitable. For the real shifts 
the analysis of (|^) has been done in [Q. For general complex shifts Ax we were unable to 
construct a Hilbert space on which both X and P would be Hermitean. We are therefore 
concentrating on two types of shifts: (a) real shifts, (b) purely imaginary shifts. 

By abuse of notation we will write Ax for real shifts, iAx for imaginary shifts, respec¬ 
tively. Thus Ax is always real. The following section is devoted to a construction of a 
suitable Hilbert space supporting imaginary shift representations. 

3. Analytic setup for imaginary shifts. We choose a particular model of a Hilbert space. 
Suppose A > 0 and consider a linear space 


Ca = {/ e L\JR), f{p) = 0 for IpI > A} 


( 8 ) 


One might think of A as a momentum cut-off and Ca is a space of wave functions compactly 
supported in the momentum space {e.g. we could fix A to be of the order of I/Ip, where Ip 
is the Plank length). In fact, we assume that A ~ 1/Ax. It is easy to see that Ca is a closed 
subspace of L^(IR). Using the Fourier transform we obtain the coordinate representation of 
Ca which will be denoted by Ha- Following |]^ we adopt the following 


Definition 1 An entire function F{x) is said to be of exponential type at most A if 

lim «nr. A 

iimsup|^l^g^ 1^1 ^ A . 

Now, one can give the following description of Ha ( Theorems 16 and 17) 


3 

























Theorem 1 (Paley—Wiener) 

H\ = {F G L^(IR) : F is entire of exponential type at most A} 

The shift operator acts on and, consequently, it acts on Ca- Using the convention 


F{x) 


1 /•+00 

_ / 

J —oo 


f{p) dp , 


we get that 


Ea, ■ f{p) ^ , / e Ca . 


( 9 ) 


( 10 ) 


Since / vanishes outside of a compact set, is unitary if = Ax G IR, Hermitean if 
Az = iAx respectively. The former case is analyzed in [Q by a slightly different method. In 
this paper we are interested in the latter case. acts on Ca as a multiplication operator 

by the function 


Dn{p) 


N 


iAx 




k=-N 


( 11 ) 


where 2 ; = For real ck^, is skew-Hermitean. In particular, this is true for the 

optimal discretization. 

Now, we turn to the study of the operator of multiplication by x. We hrst introduce a 
dense subspace of Ca defined as 

= {/ G Ca : f is a. c. , /(A) = /(-A) = o} . 

The corresponding subspace of Ha will be denoted by We note that C\ is invariant 
under E/^z- On we can integrate by parts obtaining 

dp . 

V27rJ-A I dp 


We are interested in hnding pairs X,Dn of operators satisfying Heisenberg CR in the sense 
of definition given in Sec. VI of |T^ . We look for X of the form 




— 3 z'"' 

m ^ 1 ' 1 Hm ^ 

dp dp 


( 12 ) 


Then the Heisenberg CR (|^ reads [DAf{p),X] = 1 and it implies 

\ d 




dz 


(Ba,(p)) = 1 . 


(13) 


Example For the optimal discretization scheme and N 
ai = 1/2 = CK-i. Then (^) implies 


1 we have found in 14 that 


Di{p) 


—^sinh(pAx) , 

VL^Jb 


-A < p < A , 
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and (|T^) yields 


In snmmary 


T.l^rnZ^ = 


(coshpAo:) ’ 


-A < p < A . 


-Di(p) = —^sinh(pAx) , 


X{p) = 


2i 


d d 


cosh(pAx) dp (ipcosh(pAx) 


is a conjngate pair in the sense of 

Below we will show that the pair X{p), -Di(p) is nnitary eqnivalent to the canonical 


conjngate pair 


Id 1 


i dy ’ 


tIj on 


id 


— sinhAAtc sinhAAa: 


This, in particular, implies that the 
-dimensional lattice. It is perhaps worth 


Ax ’ Ax 

spectrum of a self-adjoint extension of X{p) is a 
mentioning that we do not have in this case a doubling of canonical conjugate pairs, a 
phenomenon occurring for real shifts [|^. Finally, we would like to point out that the 
above pair X{p), Di{p) appears in |jT^. We believe that the present work provides a proper 
analytic setup for that work. For example, to dehne the position operator, the authors of 
0 had to conhne their attention to the zero momentum sector of the Hilbert space (formula 


(7) in [|T^). There is no need for such an assumption in our approach. 


4- Optimal discretization for imaginary shifts. In this Section we consider the case of 
the optimal discretization scheme. We limit ourselves to the most important aspects as the 
details are analogous to the real case studied in Sec. VI of 0- Our goal is to understand 
(P!B|), in particular, the structure of zeros of ■^Dn{p) = D'j^. We immediately have 

Lemma 1 Let Djsiip) he optimal and N he odd. Then D'^{p) has no roots forp G [—A, A]. 


Proof: The proof is very similar to the proof of Lemma 2 in |^. We compute iXxz d/dzD^ 
to obtain 


d 


(N'.r 


= (- 1 ) 

Thus zq is a zero of D'j^ iff it satishes 

(1 - z)"^^ 


N 




N I 2V \ M 
N 


vN 




z = 


-A < p < A . 


It is clear that for odd N this equation has no real roots (z > 0). This completes the proof. 
Lemma 2 If N is even and A = 1/Ax, then D'^{p) has no roots forp G [—A, A]. 

Proof: Now, the equation determining the zeros of H/f reads 


(1-^) 


2N 


2N 

N 


z = eP^^ 


-A<p< A . 


5 




















In other words 


2 pAx 1 

sinh = ajv , aw = ^ 


N 


2N 

N 


-A<p< A 


In [|T5 we have proved that {otv} is increasing and liniAr^ooaAr = 1. Note that 


0 < sinh 


2 pAx . 2 AAx 


< sinh 


Thus we have no roots iff sinh^ < oat for every even N. 
Since {otv} is increasing the latter inequality is equivalent to 


. , 2 AAx \/6 

sinh - <02 = —— 

2 4 


which implies that 


AAx < 2 arcsh 



~ 1.4379... , 


which holds if A = 1/Ax. 

Remark: It is clear from the proof of Lemma 2 that if AAx exceeds the bound 1.4379... then 
D '2 has exactly two roots. This case is then very similar to the case of real shifts ( |T^ , Sec. 
VI). From now on we assume that has no real roots on the interval of interest. One 
easily checks if —A < p < A, then iAxD'^ > 0 if iV is odd and iAxD'^ < 0 if V is even. 
Hence, the map 


a; : p — ^ ^T)w(p) , [-A, A] —> [-|T)jv(A)|, |T>7v(A)|] 


is bijective. Let us point out that iD^ip) is real for p G [—A, A]. We call the new variable 
y = iDjq{p). The corresponding Hilbert space is ([— |ZlAr(A)|, |i7jv(A)|]). The map oj 
induces an isomorphism of L^([— A,A]) and ([—|i7Ar(— A)|, |Zljv(A)|]). In particular, uj 
maps 




1 d 
i dy ’ 



(14) 


This situation should be contrasted with the results we have obtained in |]T^ for real shifts, 
for which the pairs X,Dn are, after a suitable rescaling, unitarily equivalent to two copies 
of the canonical conjugate pairs 7 ^, 7 P dehned on L^([— 1 , 1 ]). 

5. Imaginary vs. real shifts. The aim of this section is to clarify the differences between 
the real and imaginary shift representations of the Heisenberg CR. In doing so we will come 
across an interesting connection to problems in signal theory, centered around the so-called 
Shannon theorem [K|. Let us consider / G C^. Every such function can be written in terms 
of its Fourier coefficients. We dehne 


fin) = 


pA 


2A J-A 


e P J(p) dp 


(15) 
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In view of (^) we obtain 


[aj 

Since f{p) = /(’^) we conclude that 


/(") 


(16) 




Using the dehnition of F{x) we arrive at Shannon’s theorem 

FeH^ 


N sr^^/n7r\ sinfnTT — Ax) 

= ?^(t) L-a. 


(17) 


(18) 


Thus F{x) can be recovered from its sample values on the lattice 

Our objective now is to understand how the real and imaginary shifts act on those sample 
values. It is clear that the real shift for Ax = vr/A, gives (E^^F) (nvr/A) = F(^^7r). 
The question that presents itself at this point is whether there exists a simple formula for 
the action of the imaginary shift. We know that the action of Ei^x is most conveniently 
expressed on Oa by /(p) i—*• e^^*/(p). Thus, in view of ([T6|) we obtain 


and with the help of the convolution theorem and 


(eP^-/)(n) 


we get 


{EiAxF) ( 


mr\ 


\AJ 


T)"'sinh(AAx) 


■ir 


77 f mn' 


AAx 


i(n 


m)7r 


(19) 


We conclude that the imaginary shift Ei^x involves infinite number of real shifts by vr/A. 
One can easily extend these computations to the “discrete” derivative Diax- For simplicity, 
we only consider here the case N = 1 for which we obtain 


(c‘iX) (- 


'M..F) ((H) ^ (_i)»!H!^ 


i:( 


■ir 


{m — n) 71 


(AAx)2 + (m 




( m7r\ 

\~a) 


( 20 ) 


AAtc 


We remark that the main contribution to the sum (|^) comes form the region m ~ n ± 
Since A ~ 1/Ax, we conclude that (^) can be interpreted as a smeared central difference 


scheme. However, due to the slow decay one cannot ignore the tail of the expansion (pOl). It 
is perhaps fair to say that the utility of difference schemes involving inhnitely many points 
is questionable from a computational point of view. If, however, one is willing to accept that 
in principle our goal is not to hnd useful computational methods of the existing quantum 
theory but explore the foundations of quantum kinematics, then the imaginary shifts merit 
a careful study. In the next section we comment on the problem of fermion doubling viewed 
from the perspective of different choices of the discrete derivative. 

6. Remarks on fermion doubling. In this Section we briefly comment on the problem of 
the spectrum doubling for the lattice Dirac operator. For an elementary introduction to the 
problem one can consult 


12|. Let us consider the 2-D lattice Dirac equation 
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d 

^ = 75 P ^ , 


( 21 ) 


where 75 g^^Pisa momentum operator. If one uses, say, the central difference 

\E'(nAx + Ax) — \E'(nAx — Ax) 


scheme 


5 


2Ax 


followed by the plane wave substitution 

^(nAx) = , 


one obtains the dispersion relation 


u = 


sin^(pAa:) 

(Ax)^ 


Spectrum doubling arises because the right hand side has two zeros in the Brillouin zone. 
We have shown in [T^ that the spectrum doubling is tantamount to the presence of two 
copies of the canonical conjugate pairs. This effect occurs for any optimal discretization 
and real Ax. In particular, the central difference scheme is included in that scheme and 
corresponds to A = 1. The situation for the imaginary shift is different! In the simplest 
case of A = 1 we have 


P = —— sinh(pAa;) 


and the dispersion relation reads 


00 = 


sinh^(pAa;) 

Ax"^ 


-A<p< A 


-A<p< A 


There is no spectrum doubling; only a single copy of the canonical conjugate pair occurs. One 
can directly interpret (^I[) as a lattice equation by using the interpretation of the imaginary 
shift in terms of inhnitely many real shifts presented in Sec. V. 

7. Summary and conclusions. This paper concludes the project of looking for Heisenberg 
algebra representations in terms of the complex shift operators of the form (|l]) and (^. In 
the present paper we have limited our attention to the Schrodinger picture, hence to the 
Hermitean X, P pairs. For the case of purely imaginary shifts we have formulated a Hilbert 
space approach using the space of entire functions. Assuming furthermore the optimal form 


of the momentum operator introduced in |T^ we have found that (i) the spectrum of the self- 
adjoint extension of the coordinate operator X comprises a lattice, (ii) all the conjugate pairs 
are, after a suitable rescaling, unitary equivalent to the canonical pair on L^([—1,1]) 

(Hi) the pairs A, D do not exhibit the doubling which occurs for real shifts. 

As the next step within this approach one could analyze evolution (wave) equations for 
state vectors and operators corresponding to observables other than X, P. Some work in 
this direction has been done a long time ago ( [@-^). We, nevertheless, expect that the most 
interesting feature arising from our work, which is the lattice structure of the spectrum of X, 
has not been sufficiently elucidated. In particular, one could pursue a potentially interesting 
analogy with the solid state physics. But this will be a subject of the future research. 
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